High frequency dynamics in a monatomic glass. 
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The high frequency dynamics of glassy Selenium has been studied by Inelastic X-ray Scattering at 
beamline BL35XU (SPring-8). The high quality of the data allows one to pinpoint the existence of a 
dispersing acoustic mode for wavevectors (Q) of 1.5 < Q < 12.5 nm~^, helping to clarify a previous 
contradiction between experimental and numerical results. The sound velocity shows a positive 
dispersion, exceeding the hydrodynamic value by « 10% at Q < 3.5 nm~^. The dependence of 
the sound attenuation T{Q), reported for other glasses, is found to be the low-Q limit of a more 
general r{Q) oc Q{Q)'^ law which applies also to the higher Q region, where Q,{Q) cx Q no longer 
holds. 

PACS numbers: 61.43.Fs; 63.50.+x; 61.10.Eq 
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The possibility of high frequency sound propagation in 
disordered systems has stimulated a large number of in- 
vestigations in the last twenty years 0. In particular, in 
glassy systems, the existence of short wavelength acous- 
tic modes can be traced back to the pioneering numerical 
simulation works of Grest, Nagel and Rahman The 
details of these high frequency excitations are closely in- 
terwoven with some well known " anomalies" such as the 
excess of vibrational density of state (the so-called Boson 
Peak) or the thermal conductivity plateau observed in 
the 1 < r < 10 K region, which mark non-trivial distinc- 
tions from the crystalline state 0. In this respect, glassy 
selenium constitutes a system of paramount importance. 
Not only does it exhibit all the above mentioned anoma- 
lies 0, , but it is also interesting for its technological 
impact due to the well known high sensitiveness to photo- 
induced effects 0. 

Crystalline selenium exists at least in three different 
atomic arrangements. The more stable is the hexagonal 
form in which the atoms are distributed along parallel 
chains with the threefold spiral along the c axis. The 
monoclinic phases (a and b) are metastable and are both 
characterized by by eight-membered Seg-rings. Rhom- 
bohedral and cubic Se have also been reported. As the 
energy difference between the two main crystal phases 
- hexagonal and monoclinic - is very small compared 
to ksTg (Tg=303 K being the glass transition temper- 
ature), the competition between the two local structures 
endows to selenium an excellent glass-forming ability; an 
extremely rare feature for a monatomic system. For the 
above reasons, glassy selenium (g-Se) has been widely in- 
vestigated in the past by a large variety of experimental 
techniques, and also several numerical investigations have 



been reported lEHBS^I- This notwithstanding, some 
crucial aspects related to the features of the microscopic 
collective dynamics at wavelength comprising few atomic 
units are still unknown. More specifically, numerical in- 
vestigations report in this momentum transfer Q region 
evidence of collective modes up to frequencies as high as 
8 — 10 meV, characterized by a well defined phonon-like 
dispersion curve Q . However, based on recent inelas- 
tic neutron scattering (INS) experiments, a localization 
threshold of 3 — 4 mcV (Q « 3 nm^^) was proposed, 
above which the longitudinal acoustic dynamics became 
strongly overdamped [lo|. 

In this letter we present the first direct determination 
of the coherent dynamic structure factor by means of In- 
elastic X-ray Scattering (IXS) in g-Se. The extremely 
favorable inelastic/elastic intensity ratio allows for a reli- 
able analysis of the collective dynamics. We find evidence 
for an acoustic-like propagating mode with a sound ve- 
locity exceeding the hydrodynamic value, a signature of 
a relaxation process active over the probed timescale. A 
quadratic dependence of the sound attenuation on the 
mode frequency is observed all over the explored region, 
generalizing the law observed at small Q in other 
glasses. 

The experiment was performed at the high resolution 
inelastic scattering beamline (BL35XU) [ll| of SPring- 
8 in Hyogo prefecture, Japan. High resolution was ob- 
tained using the (11 11 11) reflection of perfect silicon 
crystals while a backscattering geometry (| — 6*3 « 0.3 
mrad) was used (for both monochromator and analyzers) 
in order to obtain large angular acceptance. A grazing in- 
cidence geometry for the backscattering monochromator 
was crucial to avoid heat-load broadening of the resolu- 
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tion due to the x-ray power load. The flux onto the sam- 
ple was « 3 X 10^ photons/sec (100 mA electron beam 
current) in a 0.8 meV bandwidth at 21.747 keV and a 
spot size of 80 /im diameter FWHM (full width at half 
maximum). The use of 4 analyzers crystals, placed with 
0.78 degree spacing on the 10 m two-theta arm (horizon- 
tal scattering plane), and 4 independent detectors, al- 
lowed collection of 4 momentum transfers simultaneously. 
Slits in front of the analyzer crystals limited their accep- 
tance to 0.24 nm"""^ in the scattering plane. The over-all 
resolution of the spectrometer was then 1.5 to 1.8 meV, 
depending on the analyzer crystal. Typical data collec- 
tion times were 200 s/bin, where the bin size was fixed 
at 0.25 meV. 

Elemental Selenium (purity 99.999 %) was purchased 
from Cerac Inc. The 3 mm Se shots were placed on a tem- 
perature controlled flat surface which was heated to few 
degrees above the melting point of the crystal {Tm = 494 
K). Final thick- film samples were obtained by applying 
to the viscous melt a mechanical press and fast cooling 
far below the glass transition temperature. The thickness 
of the Se specimens controlled by the pressure strength, 
lies between 50-100 /im, which is very close to the opti- 
mum for matching the absorption length of the X-rays 
(/i~^ = 50 /im). The obtained flakes where allowed to 
anneal near Tg and the IXS experiment was performed 
at T = 295 K. 

In FigQlwe report the measured spectra at selected val- 
ues of (fixed) momentum transfer. Clear evidence of an 
inelastic mode, at the wings of the elastic peak, can be ob- 
served. The behavior of the energy shift of this mode vs. 
the momentum transfer strongly resembles a phonon-like 
propagating mode. The high inelastic/elastic ratio allows 
for a detailed analysis of this mode less ar nbig uously than 
in previously studied glasses [H H H El H El [11. 
Following the prescription of generalized hydrodynamics 
E3 ) Eind modeling the vibrational dynamics as a Marko- 
vian process with an instantaneous second order memory 
function, one ends up with a damped harmonic oscilla- 
tor line shape to represent the dynamic structure factor 
S{Q,u;y. 
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FIG. 1: IXS spectra of glassy selenium (open dots) measured 
for several momentum transfers, as indicated by the labels 
(nm~^). Also reported are the instrument resolution (dot- 
ted line) and the best fit lineshape according to the model 
discussed in the text (full line). 
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The parameters il{Q) and r{Q) represent the charac- 
teristic frequency and attenuation of the mode, S{Q) is 
the static structure factor, while f{Q), the non-ergodicity 
factor, is related to the inelastic/elastic ratio. 

The measured intensity is proportional to the convo- 
lution between the instrumental resolution R{w) and the 
dynamic structure factor modified to account for the de- 
tailed balance condition, namely: 



with e{Q) being a factor accounting for the efficiency of 
the analyzers, sample thickness, the atomic form factors 
and, other angular-dependent instrumental correction fac- 
lOrs jl^. The best fitted lineshapes, obtained with Eq.|3 
are reported in Fig. ^ (full lines). It is evident that our 
model function, representing a single excitation model, 
satisfactorily accounts for the measured response. 

The results of the present IXS experiment are summa- 
rized in Fig. 121 where we report the Q-dependence of the 
three shape parameters together with that of the appar- 
ent sound velocity c{Q)=il{Q)/Q. Figure [2Ji) presents 
compelling evidence for the existence of an acoustic- 
like longitudinal branch, which is well defined up to the 
higher investigated momentum transfer (Q = 12.5nm~^). 
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FIG. 2: Relevant parameters obtained by fitting the IXS mea- 
surements, (a) Maximum of the longitudinal current cor- 
relation function (open dots) and the hydrodynamic disper- 
sion) Pi (dotted line), (b) Apparent sound velocity, Q,{Q)/Q 
(open dots), molecular dynamics results |^ (full dots), and 
the hydrodynamic value (dotted line), (c) Sound attenuation 
(open dots) with the best fit: F = O.ISQ^ with T and 
Q expressed in meV and nm~^, respectively (dotted line), 
(d) Non-ergodicity parameter as obtained from the inelas- 
tic/elastic ratio 



No evidence of localization is found at Q=3 nm~^ as it 
has been suggested in Ref.[10l|. Moreover, the disper- 
sion curve extends well beyond the boson peak frequency 
(Ebp « 1.7 meV in g-Se at T = 300 K Q), confirming 
earlier observations reported in other glasses [20j . 

The data from the INS experiment ^3 have been ana- 
lyzed assuming the simultaneous presence, in the same Q 
region under investigation here, of a longitudinal acous- 
tic (LA) mode and a local optic-like mode. The model 
that was utilized to describe the LA mode implied the 
presence of a crossover between a propagating and lo- 
calized regime, and this crossover was found to be lo- 
cated around Q « 3 nm~^ In contrast, MD sim- 
ulations have shown the existence of a well defined dis- 



FIG. 3: Comparison between the IXS results (open dots and 
error bar) and previous INS measurements from Ref. [To|. 
Also reported is the best fit to the INS data (full line) 



persion curve up to the edge of the first pseudo-Brillouin 
zone {Q « 10nm~^). This is shown in Fig[2h, where 
the sound velocity, c{Q), derived by MD (full dots) and 
by the present work (open dots) are reported. In Fig. 
13 we plot the data from the INS experiment (^^11 
dots) alongside with our IXS data (open dots) with the 
latter scaled to optimize the overlap in the common fre- 
quency range. As can be observed, the two sets of data 
are in reasonable agreement. The small differences at 
the base of the elastic peak are expected given the dif- 
ference in INS and IXS resolution functions: while both 
have Ri 1.5 meV FWHM, the IXS resolution is approxi- 
mately Lorentzian, while the INS resolution is Gaussian. 
The INS measurements, even pushing the limit of the 
technique at its maximum, span only a limited range in 
energy due to kinematic restrictions. This leads to best- 
fit lineshapes which exaggerate the 'humps' between zero 
to five meV (full line in Fig. Oj) and ultimately do not 
reproduce the IXS spectra, measured over a much wider 
energy region {—25<E<25 meV). 
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The extrapolated limit of the apparent sound velocity 
liniQ^o c(Q) « 2000 m/s (Fig.EjD), as measured by IXS, 
falls above the hydrodynamic value (cq — 1800 m/s), sug- 
gesting the presence of a possible mild positive dispersion 
effect. As recently reported for another system |21|, this 
effect may be due to a residual relaxation process -active 
in the glass- and with relaxation time in the ps range. 
The origin of this process is not clear at present, it could 
tentatively be associated with the presence of topological 
disorder which modifies the hydrodynamic sound prop- 
erties at wavelength comparable with the mean interpar- 
ticle distances l2a. 
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FIG. 4: Sound damping vs. square of the excitation fre- 
quency. The quadratic dependence on the wavevector 
shown in Fig. |2; turns out to be the low Q limit of the 
more general dependence shown here. 

The sound attenuation parameter T{Q) (Fig. Et) is 
found to be compatible with a law, up to Q's value 
where the dispersion relation is linear (i. e. c{Q) ~ con- 
stant) (Q < 5nm~^). This coincidence suggests a com- 
mon breakdown of the linear and quadratic laws for fl{Q) 
and T{Q) respectively which is expected to be caused by 
structural effects, i.e. by the Q— dependence of the static 
structure factor [l9j. In order to remove this structural 
dependence we report in Fig.^the value of sound attenu- 
ation parameter plotted against the excitation frequency 
for all measured Q values. As can be observed a simple 
power law relation now extends up to the whole explored 
momentum region. The best fitted law turns out to be 
r oc with q;=2.15±0.10. The present observation 
spotlights another aspect which needs to be encompassed 
in the explanation of the ubiquitous 
(or il{Q)^) dependence which is not yet fully explained. 

Finally we emphasize how the combination of IXS tech- 
nique and generalized hydrodynamics provides a power- 
ful tool to extract information about the non-ergodicity 



factor /(Q), i. e. the long time plateau of the inter- 
mediate scattering function once the structural arrest is 
attained. This quantity is reported Fig. Eli. It should 
be stressed that f{Q) is a true "shape" parameter (i. e. 
it does not depend on the absolute intensity, a quantity 
more difficult to access experimentally) and, in particu- 
lar, it is related to the elastic-to-inelastic intensity ratio. 
g-Se is found to have a remarkably low f{Q) value when 
compared to other glasses at Tg, a feature which allow a 
reliable lineshape analysis due to the extremely favorable 
inelastic/elastic signal. This low f{Q) value of selenium 
mirrors the high frag ility, according to the Angell clas- 
sification scheme \2^. in agreement with recent findings 
that link the value of the fragility with that of f{Q 0) 
at Tg [13. 

In conclusion, we presented an accurate measurement 
of the dynamic structure factor in g-Se, taking advantage 
of the state of the art capabilities of a new IXS facil- 
ity. The inelastic/elastic signal is particularly favorable 
in this system. Therefore, exploiting the lack of kine- 
matic constrains and the very good statistics, we have 
been able to perform a robust lineshape analysis in an 
extended and previously unexplored momentum-energy 
region. We found evidence for a well defined longitudi- 
nal acoustic mode extending all the way beyond the first 
pseudo-Brillouin zone, characterized by dispersion and 
attenuation properties typical of other glass forming ma- 
terials. In particular, the boson peak frequency lies in the 
linear dispersion region, thus suggesting a link between 
this universal feature of glasses and the reported high 
frequency acoustic excitation. Moreover, the presence of 
such a well defined mode allows to get new insight into 
the high frequency sound attenuation issue. Specifically, 
a quadratic dependence of the sound attenuation on the 
mode frequency extends the law holding in a smaller 
range and observed in other glass formers. Such an ob- 
servation indicates the prominent role of the excitation 
frequency rather than of its wavevector, and certainly 
calls for further investigations. 

The synchrotron radiation experiment was performed 
at the SPring-8 with the approval of the Japan Syn- 
chrotron Radiation Research Institute (JASRI) (Pro- 
posal No. 2003A0357-ND3-279). 
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